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Abstract 

The 0(4) supersymmetry of the hydrogen atom is utilized to construct a complete basis us- 
ing only the bound state wave functions. For a large class of perturbations, an expansion of the 
electron (exciton) wave function into such a complete set reduces the perturbed Schrodinger equa- 
tion to a standard eigenvalue problem with an equidistant unperturbed spectrum. A high-order 
polarizability of the ground state of the hydrogen atom in a static electric field is calculated via 
Rayleigh-Schrodinger perturbation theory for illustration of the new method. 
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I. INTRODUCTION 



The task of the perturbation theory is to calculate, at least approximately, the energy 
spectrum and the wave functions of the given problem using those of the already solved 
primary problem. As a rule, the spectrum of the primary problem consists of discrete and 
continuum states. Both of them should be taken into account in the expansion of the 
perturbed wave function. However, in practice the inclusion of the continuum is often an 
extremely complicated task. This is the case, for instance, of the hydrogen problem, where 
the bound state wave functions do not constitute a complete set even approximately: Some 
of the scattering states are energetically close to the bound ones, since the levels thickening 
point of the latter meets the continuum onset. 

In the present paper we show explicitly that a specific inhomogeneous coordinate scal- 
ing applied to the hydrogen bound-state wave functions converts them into a complete set. 
Early in 1935 it was shown by V.A. Fock that due to the hidden symmetry of the Coulomb 
potential (which leads to the additional spectrum degeneracy) the three-dimensional (3D) 
hydrogen problem in the momentum space can be one-to-one mapped onto a sphere or onto 
a two-sheeted hyperboloid in 4D space, for bound or scattering states, respectively- The full 
space of the hydrogen eigenfunctions is thus divided into two subspaces, the eigenfunctions 
of each subspace being an irreducible representation of the corresponding symmetry group. 
In particular, the bound-state wave functions are transformed into the hyperspherical har- 
monics, the irreducible representation of the rotational symmetry group 0(4) in the 4D 
energy-momentum space.- 

An expansion into the complete basis of hyperspherical harmonics (also referred to as 
Sturmian functions - in the coordinate representation) has been already used in nuclear 
physics when treating a three-body problem (see, e.g., Refs|3|j4j) and recently in anisotropic 
exciton problem,- when the perturbation theory was constructed on a sphere in 4D space. 
In the present work we formulate the perturbation theory directly in the coordinate space, 
basing on the complete set of functions orthogonal with the weight 1/r. The important 
feature of this basis is that only normalizable (square integrable) wave functions, which 
belong to the bound states of the perturbed problem, can be expanded. Due to the energy- 
dependent transformation of the standard Coulomb wave functions which makes the basis 
complete, the Schrodinger equation is reduced to the eigenvalue problem with equidistant 
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unperturbed energy spectrum. Finally, the Brillouin-Wigner perturbation theory converges 
for any perturbation which only renormalizes the Coulomb spectrum leaving a set of bound 
states therein. The classical example of such a perturbation is the anisotropic exciton 
Hamiltonian. In spite of its apparent simplicity and enduring interest, the exact approach 
to the anisotropic exciton problem based on the present method was developed only recently.— 
In other approaches the expansion of the wave function into an incomplete basis (see, e.g. 
Ref. |6|) led to substantial deviation in numerical results compared to the exact method, at 
least for rather strong anisotropy.- Nevertheless, in this work we do not dwell into details of 
concrete physical examples but only formulate the general recipe how to treat such kind of 
problems. 

As soon as the perturbation leads to drastic changes of the spectrum, such as decay of 
formerly bound states, Brillouin-Wigner perturbation theory is not convergent any more 
and the solution of a truncated eigenvalue problem becomes meaningless. However, finite 
orders of Rayleigh-Schrodinger perturbation theory still may have some physical meaning. 
This is the case, e.g., of the hydrogen atom (or Wannier exciton) in a static electric field 
which makes the bound states of an electron (exciton) metastable, showing up in the time 
decay of the wave functions. At the same time, the polarizability of an atom or an exciton 
is described by the second perturbation order.- To illustrate how the method works in such 
situations we calculate a high-order polarizability of the ground state in a static electric 
field. 



II. CONSTRUCTION OF THE COMPLETE BASIS 

The Schrodinger equation of the hydrogen atom has the form 

(-| v2 -7- £ ) Wr) = °- (1) 

where E, r, /i, and e are the electron energy, coordinate, mass, and the elementary charge 
(e > 0), respectively. Following Fock,- let us introduce for the bound states E < the 
energy-dependent parameter 

and make the following coordinate transformation p = r/a, where p = p ■ 
(sin # cos sin^cos^, cos9). Then, introducing a new wave function as <j>(p) oc y/pip(r), 
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we result in the following differential equation 



d d 1 
dp^ dp p 



\-He,<p) 
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where A(9, <p) is the angular part of the Laplace operator, 



+ t - A \ = 0, (3) 



X = 2^ = J h2/2 f B (4) 
a B V ~ E 

is a new dimensionless eigenvalue which originates from the Coulomb term in Eq. (CD), as = 
h 2 /pe 2 is the Bohr radius. 

Equation (j3J) has the following eigensolutions 8 

4>nlm(p) = Vnl(p) 34n(^ f), (5) 

where 3^ m are the standard spherical harmonics, {y n i} is the orthogonal set of functions 



ym( P ) = ^I^' ^-^Up) (6) 

normalized as 

oo 

Jd P y nn ( P )yM = S nnf , (7) 

n 

L q p are the generalized Laguerre polynomials (for definition see, e.g., Ref.|8|), (n, I, m) are the 
standard hydrogen quantum numbers. The eigenvalues of Eq. ([3]) have the following forrn^ 

Ail = ^ n = l, 2, (8) 

leaving the same degeneracy as of the usual hydrogen spectrum. Accordingly, the scaling 
parameter a depends on the principle quantum number as a = nas/2 that finally results in 
the standard Coulomb wave functions, normalized solutions of Eq. ([T]), 



Being the full solution of the radial wave equation following from Eq. ([3]), the set of 
functions {y n i{p}}, Eq. (jBj), also referred to as Sturmian basis,- constitutes a complete basis 
in the one- dimensional p-space. At the same time, the spherical harmonics 3^ m are the 
complete set of functions on a sphere. Consequently, the product of these two, the functions 
(f> n im(p), Eq. (jSJ), although they are not orthogonal to each other, make up a complete set in 



the three-dimensional p-space. If we now treat a as a constant parameter, one and the same 
for all (frnim, we will result in the complete set of functions in the electron coordinate space, 



X nim (r;a) = -= y nl (r/a) y lm (6, <p) = J^—if>Z m f^r) , (10) 
^ra V naa B \ la J 



normalized as 

dr 



Xt{r;a)xu'{r;a) = 6 VV >, u = (n,l,m). (11) 

Looking at the basic functions Eq. (|10p one can notice that, contrary to the Coulomb 
wave functions, all Xv have one and the same exponential decay factor e _r//2a , the remainder 
being a set of orthogonal polynomials of r. Any square-integrable function, which decays at 
r — > oo, can be expanded into the set {x^( r ; °)}- 



III. FORMULATION OF THE BRILLOUIN-WIGNER PERTURBATION THE- 
ORY 



Let us now consider the Schrodinger equation with a perturbation U which can in principle 
contain differential and/or integration operators, 



h 2 _ 9 e 2 



_^_ V 2 - - + U(r) - E 
2/i r 



tf(r) = 0. 



Expand the wave function as 



*(r) = £a(r; 



(12) 



(13) 



Let's define again the basic parameter a as in Eq. (j2J) with E being now the energy of the 
perturbed state. Then we have 



h 2 e 2 

V 2 E) Xuij-a 



h 2 1 

— (n- X)-Xnim(r;a), 



(14) 



2\i r 
where A is defined by Eq. (TjJ. 

Plugging Eq. (fl3|) into Eq. (fl2|) . using Eq. (fl4|) and the orthogonality condition Eq. (flT!) 
after convolution with Xv> w e result in the following matrix equation 



n 



n8 vv , + V vv i 



C,,i — AC,, 



(15) 



with the unperturbed spectrum defined by Eq. (jSJ) and the perturbation matrix 

2\ia 



h 2 



X*(r; a)U(r)xu'(r;a)dr. 



(16) 



At first glance, Eq. (I15p looks like a generalized eigenvalue problem, since the definition of 
V uv i contains the energy- dependent parameter a, and it seems that it cannot be solved by 
means of the direct diagonalization of the effective Hamiltonian. However, it is easily seen 
that any perturbation U can be considered as an operator which is governed by a set of 
physical parameters ati, a 2 , . . . , such as external field strength, anisotropy degree, coupling 
constant, and so on: U(r; a±, a 2 , • • • )■ After the coordinate transformation r = pa it can be 
rewritten as 

^rU(r; a u a 2 ,...) = V(p; ft, fa,...), (17) 

where the operator V does not depend explicitly on a but is a function of new parameters 
A = Piifl'i Oil, 0^2, ■ ■ ■), which, in turn, depend on a. The matrix elements Eq. f lT6|) now take 
the form [see the definition of <p u , Eq. 



J p 

Treating ft (instead of a{) as fixed parameters of the Brillouin-Wigner perturbation theory, 
we finally get rid of the a-dependence in Eq. (115)) . Thus, Eq. f|T5|) becomes the ordinary 
eigenvalue problem. The dependence of the spectrum on the physical parameters a, can 
be restored by means of inverse functions aij = «j(a; (3i, (3 2 , ■ ■ ■ )• Obviously, such back 
transformation exists only if there is one-to-one relation between a's and /3's. In other 
words, a.i should be a single-valued function. This imposes some additional restrictions 
on the perturbation U. Moreover, the requirement to have bound states in the perturbed 
spectrum implies the asymptotic behavior U(r) — > at r — > 00. In the next section we 
however consider an example of potential which does not satisfy the latter condition. 

IV. EXAMPLE: HYDROGEN ATOM IN A STATIC ELECTRIC FIELD 

The interaction of an electron with a static electric field S = (0, 0, S) has the form 

U(r) = -zeS. (19) 
According to Eq. (fT7j) . the corresponding dimensionless operator is 

V(p;f3) = -(3 p 2 cos 6 (20) 

where 
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The matrix elements V uv > are calculated in Appendix [Aj They have simple analytical form 
and satisfy rigorous selection rules. Namely, only the elements with 

m' = m, l' = l±l, n' = n, n± 1, n±2 (22) 

are nonvanishing. 

As it is already mentioned, the Brillouin-Wigner perturbation theory fails when the per- 
turbation destroys the bound states spectrum. In case of the perturbation Eq. f fT9|) . the 
electric field potential provides a non-zero probability of an electron to tunnel trough the 
local Coulomb barrier. The states become quasi-bound, having nonvanishing imaginary 
part of the energy. On the other hand, the Rayleigh-Schrodinger perturbation series is 
Borel summable, having zero radius of convergence.— Although the required accuracy can 
be achieved by sufficient number of perturbation orders taken into account, the higher the 
accuracy of calculations one wants to approach the smaller should be the region of the field 
magnitude. 

Let's calculate corrections to the electron ground state energy, accounting for only the 
first few orders of the perturbation theory. As the electron ground state is non-degenerate, 
the n-th order can be calculated by means of the following recurrent formula^ 

S V 

cf = <u C { 3) = (s > 0), 



A w - AW 



r) 
r=l 



[s > 0, v > 0), (23) 



where the upper script shows the perturbation order in the calculation of the ground state 
eigenvalue and eigenvector. For simplicity, in eigenvectors and matrix elements v — denotes 
the unperturbed ground state; it is omitted in eigenvalues A^ = Xq\ 

Due to the fact that V couples any unperturbed state with only finite number of states, 
the eigenvalues and eigenvectors of each finite order can be easily calculated analytically. 
For example, up to the eighth order, the ground state eigenvalue has the form 

A = 1 - 36 (3 2 - 20 052 (3 4 - 292 973 512 (3 6 - 73 522 89 7 716 (3 8 - . . . (24) 

To return back to the energy and electric field we use Eqs. ([2]), (@J, and (12T|) that yields 
an indirect dependence E(£). The numerical plots of the ground state energy measured in 
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FIG. 1: The ground state electron energy (measured in Rydberg) as a function of the electric field 
(measured in II£o = a B /e) calculated via the perturbation theory up to s = 2, 4, 6, 8, and 10th 
order (full curves). The same calculated as a power series of £ up to n = 2, 4, 6, and 8th order 
(dashed curves). 



Rydberg versus the electric field strength measured in £ = a 2 B /e are shown in Fig. 1, for 
different perturbation orders s taken into account. Note that even the curve s = 2 is not 
purely parabolic but contains all orders of £, since for any finite order E is a transcendental 
function of £ . Nevertheless, this function can be expressed as a finite series of £, all higher 
orders being cut off. Then we result in the known series^ 



/ £ \ 2 3 555 f £ \ 4 2 512 779 ( £ 

£n 



Ry 1 2 UJ 32 V^o 



256 



13 012 777803 fS. 

8192 \£~ ] (25) 



At first glance, this expression gives worse approximation than Eq. (j24|) . since the energies are 
shifted upwards, see dashed curves in the inset (Fig. 1). On the other hand, neither of these 
two series is convergent at higher orders and it makes sense only to speak about conditional 
convergence. As seen from Fig. 1, at £ > 0.05£o ; the series Eq. (J24l) is not convergent 
even for the lowest orders s = 2 to 10: The discrepancy in the energy calculated in two 
neighboring perturbation orders does not become notably smaller with every next order. In 
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reality, the difference between Eq. (|24|) and Eq. ( 125]) simply estimates the accuracy of the 
perturbative approach, which, in turn, is of the order of the imaginary part of the complex 
eigenenergy.— At small field strengths such decay rate of the bound electron (exciton) state 
is rather small and one can get from the perturbation theory a reasonable correction (of the 
order of several percent) to the bound state energy. 



V. DISCUSSION 



The new perturbation method formulated in the present paper has a number of advan- 
tages compared to the standard perturbation scheme. First of all, it allows to exclude from 
consideration the continuum states of the unperturbed hydrogen-like problem and, for a 
large class of perturbations, to reduce the Schrodinger equation to the ordinary eigenvalue 
problem with equidistant (instead of — Ry/n 2 ) unperturbed spectrum. In particular, the 
discrete representation of the Coulomb Green's function has the form 



G E {r, r) = ^-^Xn/ m ( r '; a)Xm m {r\ a), E < 0, (26) 

nlm 

where ME) and a(E) are given by Eqs. (j4j) and (J2j), respectively. In contrast to the results 



of Ref . 
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where the Laguerre-Sturmian expansion leads to a tridiagonal form of the Green's 
operator, the Green's function Eq. ( 126]) has the standard diagonal form with the discrete 
eigenvalues A^ -* = n snading in the denominator. 

Second, the use of the energy-dependent parameter a in the expansion of any square 
integrable function leads to one and the same exponential decay of the basic functions, the 
correct asymptotic of the expanded function. At the same time, in the standard perturbative 
methods each basic function has its own asymptotic behavior at r — > oo that basically results 
in a slower convergence compared to the present method. 

Third, the matrix elements of the perturbation operator, after it is transformed to the 
dimensionless form Eq. ( TITI) , are very simple and in many cases can be calculated analytically. 

Forth, the basic functions support the spherical symmetry that can be very useful when 
a perturbation also possesses some symmetry. For instance, the perturbative approach to 
the Stark problem, Sec. IIV] is very similar to the solution of the problem in the parabolic 
coordinates.- However, introducing of the magnetic field into the Hamiltonian, makes the 
'parabolic approach' extremely hard, since the basic functions are not the eigenfunction 
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of the angular momentum operators. Contrarily, in our method, the matrix elements of 
the electron potential in both electric and magnetic fields have simple analytical form and 
selection rules. 

Fifth, the complete set of functions Eq. ffTUl) can be used for expansion of the wave 
function even if the Hamiltonian does not contain Coulomb potential. One way to deal 
with such problems is to add and subtract the Coulomb potential from the Hamiltonian 
and and then follow the present perturbation method. Another way is to expand the wave 
function of interest into the basis Eq. (flOj) with constant (energy-independent) parameter a 
and then diagonalize the Hamiltonian.- In the latter case it is recommended to treat a as 
an adjustable parameter. 
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Appendix A: Calculation of the matrix elements 

The matrix elements of the perturbation V, Eq. (120 p . have the form 

V vv , = -(3 [ ^l(p)p 2 cos 9 MP) = -P8mm>Ju>l!L>, (Al) 

J p 

where (m = 0) 

4 = J dny lo (e)cos9y vo (9) = VW + 1 )( 2l '± V j 1 dttp^Mt), (A2) 

OO I oo 

' 

(A3) 
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Using the recurrent relations and orthogonality of the Legendre and Laguerre polynomials, 
we immediately find 

" v/(2/ + l)(2/ + 3) + * V(2Z + 1)(2Z-1) " V ' 

IXZ 1 = -V(n - l)(n + l + l)(n + l + 2)(n + l + 3) 5 n+2 , n , 
+2(2n - l)y/(n + l + l)(n + l + 2) 5 n+1 ^ 



and 



Sn^/n 2 - (I + l) 2 <5 n>n / 



+2(2n + l)y/(n - I - l)(n - I - 2) 5 n _ 1>n / 
V(w + 0(n - / - l)(n - / - 2)(n - / - 3) 6 n - 2 , n >, (A5) 



^ = (A6) 
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